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A  THERMAL  ANALYSIS  OF  HIGH-DRIVE 
RING  TRANSDUCER  ELEMENTS* 


1.  INTRODUCTION 


Heat  transfer  characteristics  are  an  important  consideration  for  designers  of  new 
underwater,  high-power  sonar  transducers.  Of  particular  concern  is  the  possible  overheating  of 
the  active  elements,  such  as  the  piezoelectric  ceramics  and  the  electrostrictive  and 
magnetostrictive  materials — all  of  which,  when  driven  at  full  power,  generate  considerable  heat. 
If  the  transducer  configuration  does  not  allow  for  adequate  cooling,  the  active  element  may  not 
perform  properly  or  reliably,  and  its  lifetime  may  be  considerably  shortened;  in  the  worst-case 
scenario,  the  active  element  may  be  destroyed.  Generally,  for  piezoelectric  ceramics,  the 
temperature  should  not  exceed  the  Curie  point  (-300  °C);  in  fact,  for  reliable  performance,  the 
temperature  should  be  held  to  levels  50%  below  the  Curie  point  (reference  2). 

Thermal  issues  are  most  often  prevalent  when  the  active  elements  are  encased  in  an 
insulating  blanket,  usually  polymer  or  a  polymer-based  composite  material.  The  fact  that  these 
materials  are  good  electrical  insulators  suggests  that  they  are  also  good  heat  insulators. 

The  thermal  performance  of  a  particular  transducer  design  is  usually  estimated  by 
performing  bench-top  testing  on  individual  components  or  by  testing  complete  units  under 
realistic  conditions.  Bench-top  testing  may  yield  data  that  do  not  correlate  with  the  thermal 
behavior  of  a  complete  unit.  Full-scale  tests,  on  the  other  hand,  should  provide  reliable  data; 
however,  this  approach  may  be  very  expensive,  especially  if  redesigns  are  required. 

This  report  considers  only  piezoelectric  ring  transducers  driven  at  high  power  levels 
under  continuous  sine-wave  drive,  which  can  create  destructive  temperatures  in  the  active 
elements.  The  transducers  consist  of  several  piezoelectric  ceramic  cylinders  of  three  sizes  that 
are  coated  with  a  polymeric  insulator  as  a  water  barrier  (see  references  3,  4,  and  5)  for 
designated  operations.  In-air  testing  of  uncoated  and  coated  ceramic  cylinders  was  first 
performed  to  determine  the  temperature-versus-time  profile  for  different  drive  levels.  The 
results  were  used  to  help  in  developing  analytical  models  that  can  predict  temperature  profiles  of 
the  transducers  in  water  and  to  help  in  providing  guidelines  concerning  thermal  issues.  The 
models  complement  the  measurements  and  may  be  used  to  optimize  material  choices  and  to 
predict  behavior  under  operational  conditions. 
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Portions  of  this  report  were  presented  at  the  137lh  Meeting  of  the  Acoustical  Society  of  America,  Berlin,  Germany, 
19  March  1999  (see  reference  1). 


This  report  focuses  on  thermal  models  for  a  compact,  high-drive,  3 1  -mode,  free-flooded 
ring  underwater  sonar  transducer.  For  this  type  of  transducer,  it  was  found  that,  thermal  issues— 
rather  than  electric  field  limits  or  mechanical  stress  limits — were  the  limiting  design  concern. 
The  models  that  were  developed  are  therefore  based  on  dielectric  losses  in  the  ceramic  material. 
The  results  are  compared  to  measurement  data  and  finite  element  (FE)  model  solutions.  These 
models  must  address  the  transient  temperature  behavior  because  the  application  requires 
knowledge  of  the  temperature  as  a  function  of  time.  Figure  1  shows  the  basic  configuration  of  a 
radially  poled,  31 -mode  piezoelectric  cylinder  that  is  addressed  in  this  report. 


Electrode 

surfaces 


R 


Figure  1.  Radially  Poled  Piezoelectric  Ceramic  Cylinder 
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2.  INPUT  ELECTRICAL  POWER 


The  radiated  acoustic  power  Pa  for  a  sonar  transducer  is  the  electrical  input  power  Pin 
minus  its  dissipated  inertial  power  losses  P/„ss;  that  is. 


P  =  P  -  P 


loss  * 


(1) 


The  electroacoustic  efficiency  is  given  by 


n  =  — — . 

lea  p 


(2) 


The  total  dissipated  power  loss  in  the  transducer  is  represented  by  the  internal  electrical  losses  Pe 
and  internal  mechanical  losses  P,„  and  is  given  by 

P^=P,  +  P.-  O) 


Total  electrical  power  input  P,„  is  thus  equal  to  the  internal  electrical  losses  Pe,  internal 
mechanical  losses  Pm,  and  radiated  acoustic  power  Pa;  that  is, 

P  =P  +  P  +P  .  (4) 

in  emu  v/ 


The  mechanical  losses  may  be  associated  with  viscous  losses,  damping  losses,  and  the 
elastic  constants  of  the  piezoelectric  material  tan  8m  (reference  6).  For  a  piezoelectric  ceramic 
cylinder,  the  viscous  and  damping  losses  are  negligible  and  the  mechanical  losses  are  associated 
with  the  elastic  constants,  which  are  typically  0.0012  for  PZT-8  and  0.0032  for  PZT-4  at  high 
electric-field  drive  conditions  (reference  6).  The  tan<5,»>  reciprocal  value  is  commonly  listed  in 
the  technical  piezoelectric  ceramics  catalogs  (as  shown  in  references  7-10)  for  piezoelectric 
properties  such  as  Qm  (mechanical  quality  factor)  and  is  typically  1000  for  PZT-8  and  500  for 
PZT-4  for  low-signal  drive  conditions  and  833  and  313,  respectively,  for  high-signal  drive 
conditions. 

A  piezoelectric  sonar  transducer  can  be  represented  by  an  equivalent  electrical  circuit  to 
describe  its  behavior  (see  figure  2).  Such  a  circuit  comprises  three  sections:  electrical 
impedance,  mechanical  impedance,  and  acoustic  impedance.  The  electrical  impedance  section  is 
composed  of  electrical  leads,  electrical  dissipative  losses  Rb ,  and  clamped  or  blocked  capacitance 
Cb  of  the  piezoelectric  ceramic  material. 

The  mechanical  impedance  section  is  composed  of  the  mechanical  vibrating  components 
such  as  mass  M  and  mechanical  compliance  Cm  of  a  radially  poled,  31 -mode  ceramic  cylinder. 
The  mechanical  losses  are  denoted  as  viscous  losses,  damping  losses  Rj,  and  the  elastic  constants 
of  the  piezoelectric  material  R,„  (reference  6). 
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The  acoustic  impedance  section  is  the  load  impedance  /?/,  exerted  by  the  medium  on  the 
radiating  surface  area  of  the  transducer.  (The  radiation  impedance  is  a  complex  quantity;  for  this 
discussion,  however,  only  the  real  part  is  addressed.)  The  electromechanical  transformer  in  the 
circuit  is  given  by  the  turns  ratio  1  :N,  which  converts  the  applied  ac  rms  drive  voltage  V  at  the 
electrical  terminals  to  an  ac  force  F  at  the  input  to  the  mechanical  section  of  the  circuit.  The 
force  drives  the  mechanical  component  and  produces  a  velocity  u  to  the  acoustic  load  /?/,,  which 
radiates  acoustic  output  power  Pa  into  the  medium. 


1  :  N 


Figure  2.  Generalized  Equivalent  Circuit  Representation  of  a  Piezoelectric  Sonar  Transducer 


2.1  LOW-FREQUENCY  APPROXIMATION 

At  low  frequencies  \/coCm  »  coM ,  where  (o  =  2rrf  is  the  angular  frequency  and  RL  «  0  is 
the  radiation  acoustic  load,  the  input  electrical  power  becomes  P,„  =  Pioss  since  Pa  =  0.  The 
equivalent  circuit  reduces  to  that  shown  in  figure  3. 


I 

- ► 


Figure  3.  Equivalent  Circuit  Representation  of  a  Piezoelectric 
Sonar  Transducer  at  Low  Frequency 


The  resulting  complex  circuit  input  admittance  is  defined  as  the  ratio  of  drive 
current-to-drive  voltage: 
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where 


Gb  —  1  /Rb  and  Cf—  Cb  +  N~Cm. 

The  input  power  at  low  frequency,  written  as  the  real  part  of  the  admittance  times  the  voltage 
squared,  is 


(6) 


The  electrical  dissipation  factor  tan  Se  at  low  frequency  is  the  ratio  of  electrical  losses  Gb  across 
the  total  capacitance  of  the  equivalent  circuit: 

tan  5  = - ! - 2 - =  (7) 

a>Rb(Cb  +  N2CJ  coC  f 

The  mechanical  dissipation  factor  tan  S„,  is  associated  with  mechanical  losses  Rm  and  mechanical 
compliance  C,„  of  the  piezoelectric  material: 


therefore,  P,„  can  now  be  expressed  in  terms  of  dissipation: 


(9) 


The  effective  coupling  factor  squared  is  defined  as 


(10) 


The  power  input  is  thus  written  as 


Pin  =  V'coC,  [tan  Se  +  k]ff  tan  8m  \ . 


(ID 


Note  that  the  rms  electric  field  in  the  ceramic  is  E  =  V/t ,  and  free  capacitance  is  defined  as 
Cf  =  K1  £(>A/t,  where  KT  is  the  dielectric  constant;  therefore. 


(12) 


V2Ct  =  E2t 2  -  =  E2KTsaAt  =  E2Kt£qVc, 

and  the  input  electrical  power  at  low  frequency  is  finally  written  as 

Pin  =o)E2K Teo Vc  tan  Se+o)E2 k)ff  K 1  so Vc  tan  8m.  (13) 

In  equation  (13),  Vc  denotes  the  volume  of  piezoelectric  ceramic. 

In  summary,  at  low  frequency,  consider  the  case  where  the  transducer’s  coupling 
coefficient  is  low:  the  major  portion  of  the  dissipated  power  P/oss  will  be  in  the  dielectric  losses, 
as  in  a  capacitor.  Thus,  for  a  ceramic  cylinder  with  a  coupling  coefficient  of  0.33,  an  electrical 
dielectric  loss  tangent  factor  of  0.0065,  and  a  mechanical  losses  tangent  factor  of  0.0012 — all 
typical  for  a  3 1  -mode  ceramic  cylinder — 98%  of  the  input  power  is  consumed  by  the  dielectric 
losses  and  2%  of  the  input  is  consumed  by  the  mechanical  losses  of  the  ceramic.  If,  however,  the 
coupling  coefficient  is  much  higher  (0.8,  for  example),  89%  of  the  input  power  is  consumed  by 
the  dielectric  losses  and  1 1  %  of  the  input  is  consumed  by  the  mechanical  losses  of  the  ceramic— 
below  resonance. 


2.2  RESONANT  DEVICE 

In  the  generalized  equivalent  circuit  (figure  2),  the  parallel  combination  of  C,„  and  Rm  can 
be  combined  into  an  impedance,  as  shown  in  equation  (14): 


R.. 


1 


1 


Z_  =■ 


joC* 


j(oCm 


j(oCm  _  jcoCn 


(1  +  7  tan  8m) 


R_  + 


jaCn 


1  + 


1 


jcoCmRn 


'  j  t3n  Sm 


1  +  tan'  5„, 


(14) 


When  tan 8m  «  1,  the  parallel  impedance  becomes  a  series  impedance: 


tan  8„, 


1 


jcoCm  coCm  jcoCn 


■+R,. 


(15) 


The  mechanical  dissipation  factor  can  now  be  represented  as 


tan  8,„  -  coR.C,„, 

m  c  tf 1  7 


(16) 


or  by  the  circuit  shown  in  figure  4. 
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Figure  4.  Parallel  Impedance  Transformed  to  a  Series  Impedance 


At  resonance,  McoCm  =  coM  and  the  equivalent  circuit  in  figure  2  reduces  to  that  in 


figure  5. 
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Figure  5.  Equivalent  Circuit  at  Resonance 
The  circuit  input  admittance  at  resonance  now  becomes 

N 2 


C  =  + 


Rc  +  R:  +  Rl 


+  j°^b- 


(17) 


The  input  power  at  resonance,  written  as  the  real  part  of  the  admittance  times  the  voltage 
squared,  is,  from  equations  (7)  and  ( 1 0) 


R_  =  V 


Gh  + 


N 


Rc  +  Rj  +  R, 


=  V  coCf  tan  Se  +  V1NzQwaCm 

=  V2coCf  tan  St,  +  V2k2ffCfcoQw 
=  V2coCf[ia.nSe+k2effQw\, 


(18) 


where  the  in-water  quality  factor  is  defined  as 
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(19) 


Q-  a>rCJR,+Ri  +  RLy 

The  in-water  loaded  quality  factor  Qw  characterizes  the  sharpness  of  the  transmitting 
pressure  response  of  the  transducer  at  its  resonant  frequency  co,  when  the  transducer  is  loaded  in 
the  water  medium.  Cm  is  the  mechanical  compliance  of  the  ceramic  cylinder.  /?/.  is  the  in-water 
load  resistance  (density  times  the  sound  speed  of  the  fluid  medium  times  the  acoustic  radiating 
surface  area  of  the  ceramic  cylinder),  and  Rc  is  the  mechanical  loss  in  the  ceramic,  as  well  as 
other  mechanical  losses  /?,  (viscous  losses,  damping  losses,  etc.). 

Consider  the  case  where  keffQw  »  1 .  The  input  electrical  power  is  approximated  as 
Pm  =  V2coC fk2ffQw  and  is  dominated  by  mechanical  losses  and  acoustic  radiation.  Alternatively, 

if  keffQw  «  1 ,  the  input  power  is  approximated  by  Pjn  =  V  2coCf  tan  Se ,  which  is  dominated  by 
electrical  losses. 

Substituting  the  electric  field  and  dielectric  constant,  equation  (13)  yields  equation  (20) 
for  the  input  power  at  resonance: 

Pm  =  q)E2Kt e0Vc  tan  St  +  coE2Q„k2ffKT£0Vc.  (20) 

The  second  term  on  the  right  side  of  equation  (20)  is  the  sum  of  dissipated  mechanical  and 
radiated  acoustic  power: 

P.  +  P^vE'-QJaK’e.K-  <2') 

The  first  term  on  the  right  side  of  equation  (20)  is  the  electrical  power  losses  of  the  piezoelectric 
ceramic  material: 

Pe  =  a>E2VcKT£0  tan&.  (22) 

The  mechanical  losses  at  resonance  are  a  function  of  the  strain  in  the  piezoelectric 
material  and  are  written  as 

Pm  ~  o)r  T2sE  Vc  tan  Sm  •  (23) 

At  low  frequencies  (that  is,  under  free  condition),  the  strain  is  S  =  dE\  but  at  resonance  the  strain 
is  amplified  by  Qw.  Hence,  S  =  dEQw,  and  the  stress  becomes  T  =  S/sh  =  EQ»d/s'\  The 
mechanical  loss,  therefore,  may  be  written  in  terms  of  the  quality  factor  and  material  constants: 

pm  =  o)rE2Q2wd2sE  Vc  tan  S„,  (24) 

where  d  is  the  piezoelectric  strain  coefficient  and  s1'  is  the  short-circuit  elastic  modulus.  Finally, 
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since  d 2  =  k2ffs E KT sa,  equation  (24)  becomes 
Pm  =  0)  E2  QlklffK  t£0  Vc  tan  5m  - 


(25) 


Table  1  lists  the  power  contribution  of  the  electrical  and  mechanical  losses  for  loaded 
quality  factors  of  3  and  30,  using  equations  (20),  (22),  and  (25),  for  an  ideal  ceramic  cylinder 
with  a  coupling  coefficient  of  0.33  and  with  other  mechanical  losses  neglected,  that  is,  Rj  =  0. 
The  dielectric  and  mechanical  loss  tangent  factor  values  in  table  1  are  at  the  maximum  electric 
field  levels  for  a  Navy  Type  Ill  ceramic  cylinder  for  these  quality  factors.  The  drive  field  levels 
are  lower  for  the  high  quality  factor  because  of  stress  limitation;  if  driven  harder,  the  ceramic 
cylinder  would  go  into  tension  and  fracture.  For  a  loaded  quality  factor  of  3,  the  mechanical 
power  losses  are  1 0  times  lower  than  the  electrical  power  losses.  The  results  reverse  for  a  loaded 
quality  factor  of  30,  where  now  the  electrical  losses  are  10  times  lower  than  the  mechanical 
power  losses  (see  also  table  8  in  reference  6). 


Table  I.  Power  Contribution  of  the  Electrical  and 
Mechanical  Losses  for  Low-  and  High-Load  Quality  Factors 


qh 

Field 

(kVrms/m) 

tanSe 

tan  Sm 

Pin 

(W) 

Pc 

(W) 

Pm 

(W) 

Pa 

(W) 

3 

400 

0.01 

0.0012 

100 

3.3 

0.35 

96.4 

30 

130 

0.005 

0.005 

100 

0.15 

14.7 

85.2 

The  data  in  table  1  demonstrate  that  internal  dielectric  losses  are  more  important  for 
lower  loaded  quality  factors  and  the  mechanical  losses  contribute  more  to  higher  load  quality 
factors.  Damping  losses  and  other  mechanical  losses,  such  as  viscous  losses,  would  change  the 
results  and  may  equalize  the  contribution  between  electrical  and  mechanical  losses  for  higher 
load  quality  factors.  The  in-water  loaded  quality  factors  for  the  ceramic  cylinders  used  in  this 
study  are  3  (as  reported  in  reference  5);  thus,  for  the  study  reported  in  the  following  section,  it 
was  assumed  that  the  mechanical  losses  are  negligible  and  equation  (22)  was  used  as  the  source 
of  dissipated  heat  in  the  analytical  temperature  models. 
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3.  ANALYTICAL  MODEL 


The  energy  balance  (heat  flow)  equation  for  a  solid  material  is  (reference  1 1 ) 


?L-aV2T+-£- 
dt  pc  ’ 


(26) 


where  T=  temperature  (°C),  t  =  time  (sec),  a=  thermal  diffusivity  (m2/s),  p=  density  (kg/m3), 
C  =  the  specific  heat  (J/kg  °C),  and  Q  =  the  rate  of  internal  heat  generation  ( W/m3).  The 
diffusivity  can  be  expressed  as 


a  = 


(27) 


where  k  =  the  thermal  conductivity  (W/m  °C). 

Conduction  in  the  solid  requires  that  the  heat  flux  in  watts, 

q  =  -kAS/T ,  (28) 

be  continuous  across  any  solid-solid  interface  of  area  A.  At  solid-fluid  interfaces  where 
convection  is  important,  the  heat  flux  can  be  expressed  as 


q  =  Ah(T-TJ,  (29) 

where  h  =  convection  film  coefficient  (W/m2  °C),  T  is  the  temperature  at  the  interface,  and  Tx  is 
the  temperature  in  the  fluid  far  removed  from  the  interface.  The  rate  of  change  of  internal 
energy  per  unit  volume  (W/m3)  is  given  by 


dU  _  dT 

- =  pC  — . 

dt  dt 


(30) 


The  internal  energy  changes  whenever  there  is  a  net  influx/outflux  of  heat  and/or  when 
internal  heat  sources  are  present.  In  this  study,  the  internal  sources  are  the  result  of  dielectric 
heating  in  the  piezoelectric  material,  as  given  by  equation  (22).  Dividing  equation  (22)  by  the 
ceramic  volume  gives  the  internal  heat  source  in  watts/volume: 

Q  =  (oE2er  tan  Se,  (31) 


where 


(32) 
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The  electric  field  may  be  written  as  V/2a,  where  V  is  the  rms  voltage  across  the 
electrodes,  and  2 a  is  the  wall  thickness  of  the  ceramic  cylinder. 

Consider  the  ceramic  cylinder  shown  in  figure  1 .  The  Laplacian  in  cylindrical 
coordinates  is  defined  as 


,  82T  1  8T  1  d2T  d2T 

T  = - —  H - 1 - ; - —  + - r". 

8r~  r  dr  r~  89'  8z~ 


(33) 


For  the  case  of  axisymmetric  geometry. 


1  82T 

r2  892 


=  0. 


(34) 


The  Laplacian  may  be  accurately  modeled  as  a  planar  region  if  the  Biot  number  is  much 
less  than  1  (reference  11).  In  this  research,  the  Biot  numbers  for  all  cylinders  studied  are  less 
than  0.05;  thus,  the  Laplacian  reduces  to 


8x~  8y~ 

where  the  Biot  number  Bi  is  a  measure 


of  the  ratio  of  convection-to-conduction,  defined  as 


(35) 


Kr2~rx) 


Bi  =  nil? — in  (36) 

K 

where  kc  is  the  thermal  conductivity  of  the  ceramic. 

The  heat  flow  equation  (26)  may  be  solved  by  the  method  of  separation  of  variables.  Let 

T(x,y,t)  =  f(t)g(x,y).  (37) 

Substituting  equation  (37)  into  equation  (35)  and  taking  a  linear  combination  of  solutions  yields 
the  complete  analytical  solution  in  the  form 


T{x,y, t)  =  X  Anm  cosf  Xn  -Icosf  Xm  ^  ][l  - e~l- '  ], 
\  a  J  V  l  ) 


K- 

a) 


(38) 


where  T„,  Am,  and  k„m  are  separation  constants  to  be  determined  from  the  boundary  conditions 
and  the  initial  conditions.  It  is  convenient  to  express  the  internal  heat  generation  term  in  the 
same  fashion  by  using  equation  (26);  that  is,  let 
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Q_ 

Pc 


cos 


(39) 


r  A 

a) 


cos  A. 


y_ 

I 


\ 

/ 


These  expansions  are  valid  provided  that  the  individual  terms  are  eigenfunctions  of  a 
linear  Hermitian  operator  (reference  12).  The  boundary  conditions  that  are  used  require  that  the 
expansions  not  be  Fourier  series;  nevertheless,  the  eigenfunctions  do  form  a  complete 
orthonormal  set  in  the  interval  -a<x<a,  -l  <y<l.  Note  that  sine  functions  are  excluded 
because  of  symmetry. 

Substituting  the  analytical  solution  equation  (38)  and  heat  source  equation  (39)  into  the 
heat  flow  equation  (26)  yields 


and 


a 


(40) 

(41) 


To  evaluate  consider  the  boundary  conditions.  Atx  =  ± a,  and  at y  =  ±l,  heat  is 

transferred  to  the  fluid  by  conduction  and  convection  according  to  equations  (28)  and  (29), 
respectively.  Thus,  at  the  interface  with  the  fluid,  continuity  of  heat  flux  is  required: 


OX 

(42) 

1 

II 

Sr- 

1 

£1 

(43) 

Applying  the  analytical  solution,  equation  (38),  into  equations  (42  and  43)  requires  that 


cot  A„  =  (— )  An, 
a  h 


(44) 


and 


cot  Am  =  (-I~r)  Am-  (45) 

/  h 

These  equations  must  be  evaluated  numerically. 

Assume  that  the  rate  of  internal  heat  generation  Q  is  uniform  in  the  ceramic.  Then,  by 
using  the  principle  of  orthogonality  of  the  eigenfunctions  in  equation  (39),  in  a  manner  similar  to 
evaluating  the  coefficients  of  a  Fourier  series,  it  is  found  that 
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_ 16  Q  sin  Xn  sin  A,„ _ 

pC  [  sin  2  +  2  ]  [  sin  2  Am  +  2  ] 


(46) 


B 


By  using  equation  (46),  the  A  coefficients  from  equation  (41)  can  now  be  determined. 
Equations  (38)  through  (46)  constitute  a  complete  solution  to  the  heat  flow  problem  that  satisfies 
the  initial  condition  T(t  =  0)  =  0. 
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4.  DISCUSSION  OF  THE  SOLUTION 


The  solution  to  the  heat  flow  equation  (38)  assumes  that  the  initial  temperature  is  zero.  If 
this  is  not  the  case  of  interest,  any  value  for  the  initial  temperature  Ta  (temperature  offset)  may  be 
added  to  the  solution  to  obtain  the  correct  result.  For  example, 

T(x,y,t)  =  X A..  cosU,  i)cos(l„  £)[1  +  T„.  (47) 

a  I 

A  representative  result  is  shown  in  figure  6.  Typically,  about  five  double  summations  in 
equation  (47)  are  enough  to  ensure  good  convergence  of  the  series. 
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Figure  6.  Temperature  Response  Profile  of  a  Ceramic  Cylinder: 
Analytical  Model,  Heater,  and  Steady-State  Conditions 


4.1  EARLY  TIME  BEHAVIOR 

The  initial  rate  of  temperature  increase  is  described  by  equation  (26).  If  a  uniform  initial 
temperature  distribution  is  assumed,  then  the  diffusion  term  is  dropped  and  equation  (26) 
becomes 


57\  _  Q 

d/',=0  pC' 


(48) 
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The  early  time  behavior  of  temperature  is  expressed  as 

T  =  -0-t  +  T.  (49) 

pC 

It  is  clear  from  the  internal  heat  source,  equations  (31)  and  (49),  that  the  initial  slope  ot 
the  temperature  profile  is  proportional  to  the  frequency,  square  of  the  applied  voltage,  dielectric 
constant,  and  loss  tangent,  and  it  is  inversely  proportional  to  the  specific  heat  C  and  wall 
thickness.  This  rate  of  heat  generation  in  the  ceramic  is  that  caused  by  dielectric  losses.  This 
portion  of  this  curve  is  referred  to  as  the  “heater.”  If  the  ceramic  material  were  perfectly 
insulated,  the  temperature  would  continuously  rise  along  this  straight  line  until  catastrophic 
failure  occurred,  as  shown  in  figure  6. 


4.2  STEADY  STATE 

The  steady-state  solution  can  be  obtained  directly  from  equation  (47)  as  the  bracketed 
exponential  function  approaches  1  and  is  expressed  by 


TSS  =  X  Anm  COS(4,  -)  COS(/Lm  y)  +  T0. 


(50) 


The  form  of  this  solution,  however,  does  not  easily  reveal  the  dependence  on  the  physical 
parameters  of  the  problem.  An  approximate  expression  that  more  explicitly  displays  the  physics 
of  the  problem  can  be  derived  by  assuming  that  the  length  of  the  cylinder  is  much  larger  than  the 
wall  thickness.  In  other  words,  the  diffusion  takes  place  primarily  radially  through  the 
cylindrical  wall  of  the  ceramic. 

The  energy  balance,  equation  (26),  may  then  be  written  for  steady  state  with 
equation  (27): 

k^  +  Q  =  0.  (51) 

dx2 


Integrating  equation  (5 1 )  and  applying  the  boundary  condition  of  equation  (42)  at  the  fluid 
interface  yields 


2k  h 


At  x  =  a,  equation  (52)  becomes 


v.s  _ 


QK 


a  I 


c  V  c 


+  T0, 


(52) 


(53) 
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where  Vc  is  the  total  volume,  Ac  is  total  surface  area  (including  inner  diameter  (ID),  outer 
diameter  (OD),  and  end  surfaces),  kc  is  thermal  conductivity,  and  a  is  one-half  the  wall  thickness 
of  the  ceramic  cylinder.  This  expression  is  an  approximate  expression  of  the  steady-state 
temperature  at  the  ceramic  wall  outer  surface. 


4.3  APPROXIMATE  TIME-DEPENDENT  MODEL 


Spatially  averaging  the  energy  balance  equations  (26)  through  (30)  over  the  volume  of 
the  ceramic  yields 


dT_ 

dt 


Q  H 

pC  pCVc 


(54) 


The  solution  to  equation  (54)  yields  the  approximate  model  as  follows: 


T  =  ( T“  -T„)  (1  -  e'/T)  +  To, 


(55) 


where  the  time  constant  to  reach  steady  state  is  r ;  that  is, 


r=PCVc 


h  Ac 


(56) 


Substitution  of  equation  (53)  into  equation  (56)  yields  the  desired  expression  for  temperature  in 
terms  of  physical  parameters: 


QVf  a  1 

Ac  U  hj 


V 


-A 
\-er 

) 


+  Tq. 


(57) 


This  approximate  model  is  in  good  agreement  with  the  exact  model  of  equation  (47)  (see 
figure  7,  which  also  shows  an  FE  solution  for  the  same  problem).  The  FE  solution  was 
calculated  by  the  FLUX2D  model  shown  in  figure  8,  which  also  provides  the  model  inputs 
(reference  13).  The  agreement  between  models  is  reasonably  good. 
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Figure  7.  Comparison  of  Exact  Analytical  Approximate ,  and  Axisymmetric  FE  Model 
Solutions  for  Uncoated  PZT-8  Ceramic  Cylinder  with  Dimensions  of  OD  =  1.826  Inches , 
ID  =  1.45  Inches,  and  Length  =  0.807  Inch  (T„  =  23  °C) 
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Figure  8.  Axisymmetric  FE  (FLUX2D)  Model  for  PZT-8  Ceramic  Cylinder  2  (See  Table  2) 
with  Dimensions  of  OD  =  1.826  Inches,  ID  =  1.45  Inches ,  and  Length  =  0.807  Inch  at 
Steady  State,  Driven  at  1000  Vrms  (Inputs  to  Model:  /=  10  kHz,  h  =  17,  K1  =  /  /  00, 

T„  =  23  °C,  Dissipation  =  0.0065,  k  =  2.1,  Cc  =  420,  pc=  7600) 
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The  approximate  model  shows  that  to  minimize  heating,  the  ratio  r/Tss  should  be 
maximized;  that  is,  r  should  be  large  and  T's  should  be  small.  This  ratio  can  be  written  as 

t/T“  =  pC  /  (58) 

a  l 

which  implies  that  a  tall,  thin  cylinder  will  dissipate  the  generated  heat  more  efficiently  than  will 
a  short  thick  cylinder. 
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5.  MODELS  FOR  COATED  CYLINDERS 


A  more  practical  case  is  a  ceramic  cylinder  encased  in  a  protective  polymer  coating  and 
driven  at  high  power  in  a  free-field  fluid  environment.  The  thermal  conductivity  and  specific 
heat  of  coatings  are  usually  less  than  those  of  the  ceramic;  therefore,  the  ceramic  temperature 
will  be  nearly  uniform  throughout  its  volume  at  any  given  instant.  This  study  treats  the  ceramic 
and  the  coating  as  lumped  parameters.  A  linear  temperature  profile  in  the  polymer  is  assumed. 
FE  models  (discussed  in  section  7)  indicate  that  the  profile  is  very  linear.  The  thermal  energy 
balance  equation  for  the  coated  ceramic  may  be  written  as 


A  C  Vc 


dT  a 

dt 


=  QVc-kPAc 


(Ta-n) 

d 


(59) 


The  thickness  of  the  coating  is  d.  The  temperature  subscript  a  refers  to  the  ceramic/polymer 
coating  interface,  and  subscript  b  refers  to  the  polymer  coating/fluid  interface  (see  figure  9). 


Figure  9.  Epoxy-Coated  Ceramic  Cylinder  Layout 


The  left  side  of  equation  (59)  is  the  total  rate  of  change  of  internal  energy  in  the  ceramic 
(equation  (30)).  On  the  right  side  of  equation  (59)  is  the  sum  of  the  total  rate  of  internal  heat 
generation  plus  the  total  heat  flux  between  the  ceramic  and  coating  due  to  conduction. 

In  the  polymer  coating,  the  energy  balance  equation  is  expressed  as 
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(60) 


^VpPpCp 


djT^  +  n) 

dt 


.  ( Ta  -_Tb )  .  L(t  _T  \ 

Ap  k p  ~  A p  ”( T h  “  T »)• 

a 


The  term  on  the  left  of  equation  (60)  represents  the  total  rate  of  change  in  internal  energy  for  the 
coating.  The  average  temperature  of  the  coating  is  represented  by  the  algebraic  sum  of  the 
temperature  at  the  ceramic  surface  and  the  temperature  at  the  coating/fluid  interface.  On  the 
right  side  of  equation  (60),  the  first  term  is  the  rate  of  heat  flow  between  the  coating  and  ceramic, 
(equation  (28)),  and  the  second  term  is  the  rate  of  convective  flow  into  the  fluid  (equation  (29)). 

Differentiating  equations  (59)  and  (60)  with  respect  to  time  can  decouple  these  two 
equations,  the  details  of  which  are  presented  in  appendix  A.  Equation  (61),  which  is  the  same  as 
that  in  equation  (A- 1 0),  represents  the  thermal  energy  balance  of  the  system  in  terms  of  the 
temperature  at  the  ceramic/polymer  coating  interface: 


d 2  Ta  ^  dTa 
,  ,  +  ei  -t— 

o  t  dt 


+  El  Ta  ~  El  Qc-> 


with 


E\  -  2  Xc  +  Tp  +  r El  -  Tc  Th\  £l~  Tp  +  r h, 


and  defining 


Qe 


,  _  Q  ..  _  M  ..  .  2M, .  . 

c  =  ~  -  J  ^  1/  '  L  P  —  J  S-  1/  '  L  h  — 


2hA, 


pcCc  '  dpcCcVc  p  dPpCpVp  PpCpVp 


(61) 


(62) 


(63) 


Here,  the  piezoelectric  ceramic  material  and  polymer  coating  material  are  represented  by  the 
subscripts  c  and  p ,  respectively.  The  Fs  and  A’s  are  the  volumes  and  areas,  respectively,  of  the 
ceramic  and  coating  (only  the  area  exposed  to  the  fluid  in  the  case  of  the  coating).  These 
parameters  have  units  of  inverse  time.  Their  inverses  represent  characteristic  thermal  relaxation 
times  for  conduction  or  convection.  For  example,  l/r<  is  a  measure  of  the  relaxation  time  for 
conduction  of  heat  from  the  ceramic  to  the  coating. 
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The  solution  to  equation  (61)  may  be  expressed  in  the  form 


Ta  = 


Ts  P,  Cc 


+  A  eXl'  +  Beh' , 


(64) 


where 


 £ 2 


£3ei 


(65) 


and 


Ai,2  ■“  £l 


i - > 

1 1 1 

1-4  — 

£.  J 

(66) 


The  plus  sign  in  equation  (66)  refers  to  A\,  and  the  minus  sign  refers  to  Ai.  The  A’s,  as  was  the 
case  for  the  r  ’s,  also  have  units  of  inverse  time.  They  are  decay  constants  for  the  overall 
system;  that  is,  their  inverses  are  characteristic  relaxation  times  for  the  complete  system.  They 
determine  the  rate  at  which  the  steady-state  temperature  is  reached. 

The  temperature  at  the  polymer/fluid  interface  may  now  be  found  from  either  equation 
(59)  or  (60)  as 


Tb  =  ( - ) 

Ts  Tc  Pc  Cc 


Ah 


+  (1  +  ^ )AeXt'  +  (1  +  —  )BeX2'. 

Tc  Tc 


(67) 


Using  the  initial  condition  that  the  temperatures  are  uniform  and  the  same  in  both  the 
ceramic,  coating,  and  the  fluid,  constants  A  and  B  can  be  evaluated  by  setting  t  =  0  in 
equations  (64)  and  (67)  and  substituting  the  initial  conditions  Ta  =  7*  =  7o  to  obtain  these  results: 


A  =  r. 


B  =  -t 


Q  (  A;  + 1 ) 

Pc  Cc  (Ai  ■  Ai) 

(68) 

Q  ( Ai  + 1 ) 

(69) 

Pc  Cc  (Ai  -  Ai) 
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6.  ASYMPTOTIC  LIMITS  AND  APPROXIMATE  SOLUTIONS 
FOR  COATED  CYLINDERS 


An  examination  of  the  initial  time  behavior  and  the  steady-state  behavior  provides  a 
means  of  verifying  the  solutions  for  coated  cylinders.  Moreover,  the  solution  is  simplified,  so 
that  the  influence  of  the  physical  parameters  is  more  apparent. 


6.1  INITIAL  TIME  BEHAVIOR 

Let  t  approach  zero  in  equation  (59)  and  integrate  with  respect  to  time.  Then,  by 
applying  the  initial  conditions  that  Ta  =  7*  =  To  at  t  =  0  and  assuming  a  uniform  temperature 
distribution,  the  same  expression  is  obtained  as  before  for  the  heater  (equation  (49));  namely, 

Ta  =  (— — )  /  +  To  as  t  ->  0,  (70) 

Pc  Cc 

which  is  in  agreement  with  the  result  expected  from  assuming  a  uniform  temperature  and 
integrating  equation  (26).  If  the  ceramic  was  perfectly  insulated,  the  temperature  would  rise 
linearly  versus  time,  that  is,  the  heater  portion  of  the  temperature  curve,  as  discussed  in  section  4. 


6.2  STEADY  STATE 


By  letting  t  approach  infinity  (oo)  in  equations  (59)  or  (60),  the  steady-state  temperature 
can  be  found  by  integrating  with  respect  to  time  and  applying  the  boundary  condition  at  the 
interface.  Assuming  Ap  «  Ac ,  at  x  =  a , 


Ta  =(}  +  —)  +  To, 

h  kP  Ac 


(71) 


where  d  is  the  thickness  of  the  coating,  and  kp  is  the  thermal  conductivity  of  the  coating  material. 
If  equation  (7 1 )  is  written  as 


Ta  -  To  =  (7  +  — ) 

h  k,,  Ac 


(72) 


it  can  be  seen  that,  for  a  given  volume-to-area  ratio  and  a  given  heat  source  Q ,  the  controlling 
parameters  are  \/h  and  d/kP .  If  h  is  very  large,  as  is  the  case  for  water  where  h  is  on  the  order  of 
100  or  greater,  the  steady  state  is  determined  by  the  thickness-to-conductivity  ratio  {d/kp).  If, 

as  is  the  case  for  air,  h  is  relatively  small,  then  the  convection  process  dominates.  When  the 
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length  of  the  ceramic  is  much  greater  than  its  wall  thickness,  Vj  Ac  »  a ,  equation  (72)  can  be 
approximated  as 


Ta  -  To  =  (7  +  — )  «  Q, 
h  kn 


(73) 


which  illustrates  a  linear  dependence  of  the  steady-state  condition  on  the  thickness  of  the 
ceramic. 

6.3.  SIMPLIFIED  SOLUTIONS 

The  general  solution  for  a  coated  cylinder,  equation  (67),  contains  two  time  constants,  A\ 
and  A2.  In  most  instances,  these  differ  by  several  orders  of  magnitude.  Thus,  one  component 
vanishes  very  early  in  time  relative  to  the  longer-lived  component.  Equation  (64)  may  thus  be 
rewritten  as 


Ta 


Q 


+  Be*2'  l  +  —e(X'  h)>  +To, 

B 


(74) 


ts  Pc  Cc 


or,  substituting  equations  (68)  and  (69), 


- s -  +  Be*2'  1 

ts  Pc  Cc 


(75) 


Both  time  constants  are  negative  (as  they  should  be  physically),  but  one  can  see  from 
equation  (66)  that  U,|»U2|.  Thus,  the  solution  can  be  approximated  as 


To  =  To  + - - -  +  BeklJ. 

ts  Pc  Cc 


(76) 


Using  equation  (69)  yields 


(77) 


Recall  that  at  initial  time  t  =  0,  Ta  =  T0,  which  implies  that 


in  equation  (75)  and,  thus, 

Ta  *  - ^ -  (1  -  ehJ)  +  To-  (79) 

Ts  Pc  Cc 

Figure  10  compares  the  simplified  model  (equation  (79))  with  the  more  exact  result 
obtained  with  equation  (64).  The  temperature  profile  computed  with  a  FLUX2D  FE  model  is 
included  in  figure  10  and  is  also  shown  in  figure  1 1  (model  inputs  are  also  included  in  figure  1 1). 


Figure  10.  Comparison  of  Exact  Analytical  Solution  and  Approximate  Solution  for  Coated 
PZT-S  Ceramic  Cylinder  3  (See  Table  2)  with  Dimensions  of  OD  =  1.054  Inches, 

ID  =  0. 791  Inch,  and  Length  =  0.527  Inch  (T„  =  25°C) 
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Figure  11.  Axisymmetric  FE  (FLUX2D)  Model  for  Coated  PZT-8  Ceramic  Cylinder  3  (See 
Table  2)  with  Dimensions  ofOD  =  1.054  Inches,  ID  =  0. 791  Inch,  and  Length  =  0.527  Inch 
at  Steady  State,  Driven  at  1000  Vrms  (Input  to  Models:  /=  10  kHz,  h  =  22,  Kr  =  11 00, 

T„  =  25  °C,  Dissipation  =  0.01,  k  =  2.1,  Cc  =  420,  pc  =  7600) 
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7.  MEASUREMENT  VERSUS  MODEL 


7.1  BARE  CYLINDERS  IN  AIR 

The  first  thermal  tests  were  conducted  on  bare  cylinders  in  air.  The  temperature  was 
measured  by  a  thermocouple  attached  to  the  outer  radius  of  the  cylinder  at  the  midpoint  of  the 
length  (see  figure  12).  Two  important  quantities  were  the  (1)  temperature-versus-time  profile  of 
the  outer  surface,  and  (2)  temperature  at  the  region  midway  between  the  inner  and  outer  surfaces. 
The  center  region  is  expected  to  be  the  hottest,  but  there  is  no  practical  way  to  measure  the 
temperature.  The  models  developed  in  section  4  were  used  for  comparison.  (The  models  are 
analytical  models  (that  is,  closed-form  equations),  but  some  finite  element  analyses  (FEA)  were 
also  used  to  verify  the  models  and  the  data.) 


Temperature  Probe 
(Thermocouple  wire) 


31 -mode 

Ceramic  Cylinder 


Drive  wires 


Figure  12.  Uncoated  Ceramic  31 -Mode  Cylinder  with  Thermocouple  Wire  Epoxied 


Three  cylinders  with  different  dimensions  (see  table  2)  and  properties  (see  table  3)  were 
modeled.  These  cylinders  were  suspended  in  air  by  their  drive  wires  in  an  enclosure  to  reduce 
air  currents.  Note  that  one  cylinder  was  composed  of  PZT-4  (Navy  Type  I)  material  and  the 
others  of  PZT-8  (Navy  Type  Ill)  material.  Each  cylinder  was  driven  at  three  voltage  levels  (500, 
750,  and  1000  Vrms)  below  their  first  resonance  while  their  electrical  input  current,  phase,  and 
power  were  monitored.  The  modeling  results  for  these  three  cylinders  are  shown  in  figures  13, 
14,  and  15. 
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Table  2.  Piezoelectric  Ceramic  Cylinder  Dimensions  and  Measured  Small  Signal  Properties 


Cylinder 

No. 

OD 

(in./cm) 

ID 

(in./cm) 

Length 

(in./cm) 

Material 

Resonant 
Frequency 
(kHz)  ‘ 

keff 

Cap 

w 

Diss 

(tan<5E) 

1 

2.546/ 

6.467 

2.160/ 

5.486 

0.964/ 

2.449 

PZT-4 
(Type  1) 

17.2 

0.34 

10.8 

0.0032 

2 

1.826/ 

4.638 

1.454/ 

3.693 

0.807/ 

2.050 

PZT-8 

(Type 

III) 

26.2 

0.33 

5.5 

0.0026 

3 

1 .054/ 
2.677 

0.791/ 

2.009 

0.527/ 

1.339 

PZT-8 

(Type 

_ uu _ 

47.1 

0.32 

2.5 

0.0034 

Table  3.  Piezoelectric  Ceramic  Material  Properties 
(References  7-10) 


Property 

PZT-4  (Navy  Tvpe  1) 

PZT-8  (Navy  Type  III) 

Density  p(kg/m3) 

7600 

7600 

Curie  Temperature  (°C) 

>300 

>300 

Heat  Capacity  Cc  (J/kg  °C) 

420 

420 

Thermal  Conductivity  kc  (W/m  °C) 

2.1 

2.1 

Dielectric  Constant  Kr  ( e '  /en ) 

1300 

1000-1100 

Dissipation  Factor  (tan<S) 

Large  Signal  (4  kVrms/cm) 

0.04 

0.01 

Figure  13.  Modeled  (+,  o,  *)  and  Measured  ( — )  Uncoated  Cylinder  I 

In- Air  Results  at  5  kHz 
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Figure  14.  Modeled  (+,  o,  *)  and  Measured  ( - )  Uncoated  Cylinder  2 

In- Air  Results  at  10  kHz 
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Figure  15.  Modeled  (+,  o,  *)  and  Measured  ( - )  Uncoated  Cylinder  3 

In- Air  Results  at  10  kHz 


The  modeled  rate  of  internal  heat  generation  per  unit  volume  was  assumed  equal  to  the 
full  measured  input  wattage  at  each  measured  temperature.  To  validate  this  assumption,  the 
internal  heat  generated  was  also  computed  from  equation  (31)  using  measured  values  of 
dissipation  and  dielectric  constants  (see  appendix  B).  Note  that  measured  values  were  used 
because  the  “book  values”  (references  7-10)  for  permittivity  and  dissipation  did  not  correspond 
to  measured  values.  In  each  case,  both  sources  produced  identical  time-temperature  curves. 

Note  that  the  curves  generated  by  the  model  are  not  smooth — caused  by  a  non-smooth  variation 
in  the  dissipation  curves  in  appendix  B.  This  variation  in  the  dissipation  may  be  attributed  to 
measurement  error  resulting  from  the  measured  quantities  being  near  the  resolution  limits  of  the 
instrument.  The  dielectric  constant  also  varies  with  temperature,  but  quite  smoothly.  The 
measured  electrical  input  power,  dielectric  constant,  and  dielectric  dissipation  factor  for  each 
cylinder  are  shown  in  appendix  B.  A  measured  example  of  a  cylinder  in  thermal  runaway 
leading  to  fracture  of  the  cylinder  is  shown  in  figure  16. 


0  5  10  15  20  25  30 

Time  (min) 

Figure  16.  Measured  Example  of  a  Cylinder  in  Thermal  Runaway 
Leading  to  Fracture  of  the  Cylinder 


Variations  in  the  film  coefficient  h  show  a  consistent  pattern:  they  decrease  with  cylinder 
size  and  increase  with  frequency  (see  table  4).  For  engineering  purposes,  the  film  coefficient 
depends  primarily  on  the  inverse  of  the  boundary  layer  thickness  (reference  1 1).  Thus,  it  is 
reasonable  to  see  a  variation  of  h  with  cylinder  size;  however,  the  dependence  on  frequency  is 
difficult  to  explain.  An  in-depth  analysis  of  the  film  coefficient  is  beyond  the  scope  of  this 
report,  but  the  values  of  h  used  are  in  the  range  of  values  used  in  engineering  calculations  for  the 
convection  of  air  over  smooth  surfaces  (reference  14)  (see  table  5).  Choosing  h  in  water  was 
conveyed  via  private  communication  with  a  heat  transfer  specialist  (reference  15). 
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Table  4.  In- Air  Convection  Film  Coefficients  (It)  for  Different  Ceramic  Cylinders 


Cylinder  No. 

5  kHz 

10  kHz 

15  kHz 

1 

10 

— 

— 

2 

— 

17 

22 

3 

— 

22 

30 

Table  5.  Representative  Values  of  Convective  Heat  Transfer  Coefficients  (Approximate) 


Type  of  Convection:  Fluid 

h  (W/m2  °C) 

Free,  Air 

5-25 

Forced,  Air 

10-500 

Forced,  Water 

100-  15,000 

Forced,  Boiling  Water 

2500  -  25,000 

Forced,  Condensing  Steam 

5000-100,000 

7.2  COATED  CYLINDERS  IN  AIR 

The  same  set  of  cylinders  was  coated  with  an  epoxy,  and  the  same  set  of  measurements 
were  taken  as  those  for  the  uncoated  set  of  cylinders.  The  applied  epoxy  has  a  relatively  high 
thermal  conductivity  (k»  1)  and  high  heat  capacity  (C  «  1400),  but  a  low  electrical  conductivity. 
Thus,  the  epoxy  was  expected  to  provide  adequate  water  barrier,  electrical  insulation,  and  good 
thermal  properties  for  the  removal  of  heat  from  the  ceramic  core,  which  is  vitally  important  for 
high-power  drive  in  water.  Polyurethane  and  epoxy  coating  material  properties  are  given  in 
table  6.  Figure  17  shows  a  representative  coated  cylinder  (figure  9  in  section  5  presents  the 
layout  used  to  model  this  case).  Measurements  were  compared  with  results  using  equation  (64) 
in  figures  1 8  through  20.  The  values  of  h  assumed  in  the  model  were  the  same  as  those  used  for 
the  bare  cylinders  (see  table  4).  For  each  cylinder,  at  each  drive  level  tested,  the  analytical 
model  predicts  the  measured  temperature  rise. 


Table  6.  Coating  Materia!  Properties 


Property 

Polyurethane 

Epoxy 

Density  (kg/m3) 

1100-1200 

1830 

Heat  Capacity  (J/kg  °C) 

80 

1400 

Thermal  Conductivity  (W/m  °C) 

0.2 

1 

Sound  Speed  (m/s)  Bulk/Bar 

1400-1700/NA 

2700/1900 

N/A  =  not  applicable. 
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Epoxy  Coated  31 -mode 
Ceramic  Cylinder 


Figure  17.  Epoxy-Coated,  31 -Mode  Ceramic  Cylinder  with  Thermocouple 
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Figure  18.  Modeled  (+,  o)  and  Measured  ( - )  Coated  Cylinder  1 

In-Air  Results  at  5  kHz,  Coating  Thickness  =  0.023  Inch  (0.58  mm) 
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In-Air  Results  at  10  kHz ,  Coating  Thickness  =  0.022  Inch  (0.56  mm) 
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In- Air  Results  at  10  kHz,  Coating  Thickness  =  0.021  Inch  (0.54  mm) 


35 


7.3  IN-WATER  MEASUREMENTS  COMPARED  WITH  MODELS 


A  thermal  life  test  was  conducted  on  a  free-flooded  cylinder  (cylinder  2)  coated  in  a 
0.022-inch-thick  layer  of  a  thermal  conductive  epoxy  (E-301  AD,  manufactured  by  CASTALL 
Inc.,  Weymouth,  MA).  A  50-foot  thermocouple  wire  (temperature  probe)  was  embedded  in  the 
epoxy  coating,  and  the  temperature  was  monitored  on  a  digital  temperature  meter,  as  shown  in 
figure  17.  Testing  was  conducted  at  the  Acoustic  Pressure  Tank  Facility  (APTF)  at  the  Naval 
Undersea  Warfare  Center  (NUWC)  Division,  Newport,  RI. 

The  APTF  measurement  test  vessel  is  a  closed,  water-filled  steel  tank  that  measures 
36.5  feet  long  by  12.5  feet  in  diameter  and  can  be  pressurized  to  2700  psig.  The  free-flooded 
cylinder  was  positioned  in  the  center  of  the  tank,  and  the  hydrophone  was  positioned  at  a 
4.15-meter  distance  for  monitoring  the  generated  sound  pressure  level.  The  tank  was  pressurized 
to  50  psig  to  avoid  cavitation  that  would  be  generated  by  the  transducer  at  a  high  sound  pressure 
level.  The  cylinder  was  tested  with  three  different,  constant,  steady-state  drive  voltage  levels  at 
26  kHz  for  20-minute  durations.  The  cylinder  was  allowed  to  cool  in  the  tank  to  ambient 
temperature  ( 1 6.6°C)  between  changes  in  drive  voltage.  The  drive  voltage  levels  were 
1000  Vrms,  1400  Vrms,  and  1750  Vrms,  which  equate  to  source  levels  of  193  dB,  196  dB,  and 
198  dB  re  1  pPa  at  1  meter,  respectively  (see  the  measured  transmit  voltage  response  (TVR)  of 
cylinder  2  plotted  in  figure  21).  As  shown  in  the  TVR,  a  free-flooded  cylinder  produces  a 
double-resonant  transmitting  pressure  response,  where  the  lower  resonance  is  generated  by  the 
inner  cavity  of  the  cylinder  acting  as  a  Helmholtz  resonator,  and  the  upper  resonance  is  caused 
by  the  fundamental  circumferential  cylinder  resonance  that  radiates  acoustic  energy  from  both 
the  inner  and  outer  surfaces  of  the  cylinder  walls.  This  cavity  resonance  is  combined  with  the 
“ring”  resonance  to  produce  a  broadbandwidth  transmitting  response  that  is  depth  independent. 


Figure  21.  Measured  TVR  of  Cylinder  2 
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Absolute  sound  pressure  level  measurements  could  not  be  made  in  the  tank  because  of 
reflecting  and  standing  waves  caused  by  the  steady-state  continuous  drive  signals  that  were  used. 
The  measured  electrical  input  powers  associated  with  the  three  different  drive  voltages  are  shown 
in  figure  22. 
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Figure  22.  Measured  Electrical  Input  Power  Levels  of  Free-Flooded  Test  Cylinder  for  Three 
Different  Steady-State  Voltage  Drive  Levels  at  26  kHz  (Tank  Pressure  50  psig) 


The  relative  flatness  of  the  power  responses  over  time  indicate  that  the  impedance 
magnitude  and  impedance  phase  angle  of  the  cylinder  remain  constant  over  time  and  that  heat 
energy  generated  in  the  cylinder  (caused  by  the  electrical  drive  power)  is  dissipated  through  the 
epoxy  into  the  surrounding  water  as  fast  as  it  is  generated.  This  finding  is  evident  by  the 
measured,  in-water  steady-state  temperature  profile  for  the  different  voltage  drive  levels  (see 
figure  23).  The  in-water  analytical  results  for  coated  cylinder  2  is  shown  in  figure  24  for 
1400-Vrms  drive  and  a  film  coefficient  of  1350  W/nr  °C. 
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Figure  23.  In-Water  Measured  Cylinder  2  with  Epoxy  Coating  at  26  kHz 
(Ring  Resonance)  and  1400-Vrms  Drive 


Figure  24.  In-Water  Analytical  Model  of  Cylinder  2  with  Epoxy  Coating  at  26  kHz 
(Ring  Resonance)  and  1400-Vrms  Drive 
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For  validation  of  the  simple  analytical  model  and  measurements,  a  coupled 
electromechanical  thermal  analysis  was  performed  with  the  AT1LA  FE  code.  In  the  AT1LA 
code,  losses  in  active  and  passive  materials  are  accounted  for  by  using  complex  physical 
constants  for  dielectric,  mechanical,  and  piezoelectric  properties.  The  dissipated  power  for  each 
element  and  for  the  whole  transducer  structure  can  be  obtained  (reference  1 6);  a  thermal  analysis 
by  heat  transfer  that  includes  conduction  and  convection  provides  the  temperature  profile  in  the 
transducer.  For  a  steady-state  solution,  thermal  behavior  is  weakly  coupled  to  the 
electromechanical  response.  The  solution  may  thus  be  performed  as  a  two-step  analysis,  for 
which  the  electromechanical  behavior  is  first  computed  and  the  resulting  dissipated  power  is 
then  applied  as  a  heat  generator  to  determine  the  resulting  temperature  of  the  system. 

Figure  25  compares  the  measured  and  computed  ATILA-modeled  acoustic  TVR  for 
cylinder  2.  The  modeled  acoustic  performance  closely  matches  the  measured  response. 

Figure  26  depicts  cylinder  2  with  a  thermally  conductive  epoxy  coating  that  was  driven  in  water 
at  the  ring  resonance.  The  FE  model  tracks  the  steady-state  solution  temperature.  The 
temperature  profile  computed  at  the  ring  resonance  is  depicted  in  figure  27,  where  a  section  of 
the  cylinder  and  epoxy  coating  is  hidden  to  display  the  temperature  through  the  cylinder 
thickness.  Note  the  temperature  at  equilibrium  remains  near  ambient  because  of  the  good 
volume-to-surface-area  ratio  for  the  free-flooded  projector.  If  the  temperature  rise  was 
significant,  however,  because  material  properties  are  temperature  dependent,  additional  analysis 
steps  would  be  performed  with  updated  material  properties  until  equilibrium  is  reached. 


Figure  25.  Comparison  of  Measured  and  FEA-Modeled  (A  TILA)  Acoustic  TVR 

for  Free-Flooded  Cylinder  2 
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Figure  26.  Measured  and  Steady-State  (SS)  FEA  Mode l  of  Free-Flooded  Cylinder  2  in  Water , 
Driven  at  1000  Vrms  and  1400  Vrms  at  26  kHz  (Ring  Resonance) 
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Figure  27.  In-Water  Steady-State  A  TILA  FEA  Model  Temperature  Profile  Through  a  Section 
of  Cylinder  2 ,  with  Epoxy  Coating  Cut  Back ,  Driven  at  1400  Vrms  at  26  kHz  (Ring 
Resonance ),  Coating  Thickness  d  =  0.022  Inch  (Inputs  to  Model:  f—  26  kHz , 
h  =  1350  W/m2  °C,  K7  =  1100 ,  T()  =  16.5  °C,  Dissipation  =  0.01 , 
k  =  2.1,  Cc  =  420,  pc=  7600) 
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7.4  COMPARISON  OF  COATINGS 


The  analytical  model  represented  by  equation  (6 1 )  is  based  on  the  assumptions  that  the 
temperature  in  the  ceramic  is  uniform  and  that  the  temperature  varies  linearly  in  the  coating.  To 
verify  this  assumption,  an  FE  model  was  used  to  calculate  the  temperature  profiles  through  the 
cross  section  of  the  ceramic  and  epoxy  and  polyurethane  coatings.  Figure  28  shows  the  in-air 
steady-state  solution  for  the  FE  model  calculated  with  both  the  high  thermal  conductivity  epoxy 
coating  used  in  the  measurements  and  a  lower  thermal  conductivity  polyurethane.  The  results  of 
this  FE  model  demonstrate  that  the  assumptions  are  reasonable. 


Figure  28.  In- Air  FE  Temperature  Contour  of  Cylinder  2  from  Center  of  Cylinder  x  to  Outer 
Layer  of  Coating  d  with  and  without  Coating ,  Driven  at  1000  Vrms  (a  =  Outer  Radius  of 
Cylinder,  T„  =  23  °C),  Coating  Thickness  d  =  0.022  Inch  (Inputs  to  Model:  f=  / 0  kHz,  h  =  17 
W/m2  °C,  Kt  =  1100,  T„  =  23  °C,  Dissipation  =  0.01,  k  =  2.1,  Cc  =  420,  pc  =  7600) 


Comparing  the  results  for  the  coated  cylinders  with  those  for  the  uncoated  cylinders 
reveals  that  the  temperature  profiles  are  about  the  same  for  these  thin  coatings.  This  result  is  a 
consequence  of  an  earlier  observation  (in  section  6.2)  that  when  the  l/h  term  dominates  ( h  small) 
over  the  d/kp  term,  the  heating  process  is  dominated  by  convection.  In  the  case  of  holding  the 
convection  coefficient  constant,  the  high  conductivity  coating  will  draw  heat  away  from  the 
ceramic  cylinder  (that  is,  it  will  act  as  a  heat  sink). 

For  comparison,  in-water  FE  models  were  used  to  calculate  the  steady-state  temperature 
profiles  in  the  coatings  as  the  thickness  was  increased  to  0.2  inch,  which  is  typical  of  a 
commercially  available  free-flooded  ring  transducer  coated  with  polyurethane.  Figure  29  shows 
the  results  for  an  ambient  water  temperature  of  20°C  of  cylinder  2  driven  at  1000  Vrms.  The 
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results  show  that  the  polyurethane  diverges  at  a  much  greater  slope  than  does  the  epoxy  coating 
when  compared  to  a  non-coated  case.  The  epoxy  coating  thickness  can  be  significantly  thinner 
than  the  polyurethane  coating  because  it  is  a  tougher  material.  Comparing  the  epoxy  coating 
thickness  of  0.0625  inch  with  polyurethane  coating  thickness  of  0.2  inch,  the  polyurethane 
temperature  is  30°C  higher.  If  the  water  temperature  was  40°C,  as  encountered  in  various  ocean 
environments,  and  the  drive  voltage  of  the  cylinder  was  1750  Vrms  (as  in  section  7.3),  the 
polyurethane-coated  cylinder  would  reach  a  temperature  of  138°C  (close  to  one-half  the  Curie 
temperature  of  PZT)  and  an  epoxy-coated  cylinder  would  reach  a  temperature  of  50°C  for  a 
0.0625-inch  thickness  or  63°C  for  a  0.2-inch  thickness. 


Figure  29.  In-Water  FE  Steady-State  Temperature  Profile  of  Cylinder  2  as  Thickness  of 
Epoxy  and  Polyurethane  Coatings  Is  Increased  (Inputs  to  Model:  f=  26  kHz, 
h  =  1350,  Kt  =  1100,  T„  =  20 °C,  Dissipation  =  0.01,  k=2.1,Cc  =  420,  pc  =  7600) 

Table  7  shows  the  FE  model  calculated  temperatures  for  different  coating  thicknesses  and 
different  drive  levels  for  the  epoxy  and  polyurethane  coatings.  Note  that  the  temperature 
increase  with  drive  level  is  linear  with  the  voltage  squared  (or  electric  field),  as  in  equation  (3 1 ). 
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Table  7.  FE  Model  Calculated  Steady-State  Temperature  for  Different  Coating  Thickness 
(Inputs  to  Model:  f=  26  kHz,  h  =  1350,  Kr  =  1100,  T„  =  0°C, 

Dissipation  =  0.01,  k  =  2.1,  Cc  =  420,  pc  =  7600) 


Coating 

Temperature  (°C) 

Thickness 

Epoxy 

Polyurethane 

(inch) 

1000  Vrms 

1400  Vrms 

1750  Vrms 

1000  Vrms 

1400  Vrms 

1750  Vrms 

0.0 

1.4 

2.7 

4.1 

1.4 

2.7 

4.1 

0.22 

2.0 

4.0 

5.9 

5.0 

9.9 

15.9 

0.0625 

3.3 

6.5 

10.2 

11.0 

23.7 

36.9 

0.2 

7.6 

14.9 

23.3 

32.0 

62.7 

98.0 

Figure  30  compares  the  in-water  temperature  profiles  of  the  two  coatings  different 
thermal  properties.  The  epoxy  coating  has  a  high  thermal  conductivity  and  heat  capacity;  the 
polyurethane  coating  is  a  polymer  whose  thermal  properties  are  two  orders  of  magnitude  lower 
than  those  of  the  epoxy. 


thi(  rtn  ii) 


Figure  30.  In-Water  Model  of  Cylinder  2  with  Epoxy  Coating  and  Polyurethane  Coating 
at  26  kHz  (Ring  Resonance)  and  1000-Vrnts  Drive,  Coating  Thickness  d  =  0.0625  Inch 
(Inputs  to  Model:  f=  26  kHz,  h  =  1350,  Kr  =  1  WO,  T„  =  20% 7, 

Dissipation  =  0.01,  k  =  2.1,  Cc  =  420,  pc  =  7600) 
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Recall  from  equation  (72)  that,  when  h  is  large  (as  in  water),  the  steady  state  is  essentially 
determined  by  the  thickness-to-conductivity  ratio;  that  is. 


rss 
a 


vc  Q 


Ac 


(80) 


The  relative  performance  of  the  two  coatings  (epoxy  and  polyurethane)  may  then  be  expressed  as 
the  ratio  of  their  conductivities: 


rss  rr>  I 

p  ~  1  o  _  Ke_  /g  j 

Te  -  To  kP 

Thus,  for  a  given  thickness  d,  the  steady-state  temperatures  (minus  Tt))  in  water  for  two  coatings 
of  conductivity  kE  and  kp  are  related  by  the  ratios  of  their  conductivities.  The  coating  materials 
described  in  table  6  have  a  conductivity  ratio  of  5,  with  the  steady-state  values  (minus  T(>)  in 
figure  30  approaching  the  relationship  in  equation  (81). 
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8.  CONCLUSIONS 


The  overall  goal  of  this  project  is  to  characterize  and  alleviate  overheating  issues  in  the 
active  driving  piezoelectric  ceramic  elements  and  the  associated  passive  materials  used  in 
transducer  designs  for  high-power,  compact  applications.  The  investigation  strives  to  identify 
design  criteria  that  avoid  reduced  acoustic  performance,  reliability,  and  service  life  of  a 
transducer  caused  by  thermal  limits.  Analytical  and  FE  models  provide  the  means  to  evaluate 
thermal  behavior  during  the  design  stage  of  new  projector  development. 

The  specific  objective  of  the  work  reported  herein  was  to  develop  analytical  methods  to 
evaluate  the  thermal  behavior  of  free-flooded,  3 1 -mode  cylinder  transducers.  The  principles 
developed  in  this  report  may  be  extended  to  other  sonar-type  transducers,  such  as  Tonpilz  and 
flextensional.  The  analytical  results  were  confirmed  by  both  FE  methods  and  measurements. 

The  simplicity  of  the  free-flooded,  31 -mode  cylinder  transducer,  which  has  both  a  high 
mechanical  quality  factor  in  air  and  a  low  quality  factor  in  water,  provided  a  convenient  design 
for  this  investigation  because  it  permits  evaluation  of  the  impact  of  material  losses  on  heat 
generation  both  below  resonance  and  at  resonance.  Through  measurement  and  model  prediction, 
the  temperature  rise  in  the  piezoelectric  material  properties  was  dominated  by  the  dielectric 
dissipation  below  and  at  resonance  for  the  case  of  free-flooded,  31 -mode  cylinder  transducers 
driven  in  water. 

This  report  demonstrates  how  thermal  behavior  can  be  modeled  for  some  simple 
geometries  using  one-dimensional,  partial-differential  equations.  This  analytical  model  may  not 
apply  to  all  cases  of  interest;  the  intent  is  to  demonstrate  the  general  problem  and  to  ultimately 
gain  a  better  understanding  of  thermal  behavior. 

The  steady-state  solution  of  the  FE  model  captures  the  temperature  profile  throughout  the 
projector  at  equilibrium.  A  steady-state  solution  may  be  used  to  evaluate  the  temperature 
resulting  for  a  given  duty-cycle  drive  by  reducing  dissipated  power  density  by  the  percent  duty 
cycle.  For  a  transient  solution,  the  thermal  conduction  and  heat  capacity  of  the  materials  must  be 
included  in  an  FE  system  of  equations.  The  fully  coupled  system  may  then  be  solved  for  each 
time  interval  using  a  method  of  weighted  residuals.  The  temperature  rise  in  a  transducer  also 
causes  a  thermal  stress,  which  may  be  taken  into  account  using  an  iterative  procedure,  which  is 
the  subject  of  another  investigation. 
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APPENDIX  A 

DECOUPLING  OF  THERMAL  ENERGY  BALANCE  EQUATIONS 
FOR  COATED  CYLINDERS 


The  energy  balances  for  a  coated  ceramic  cylinder  and  the  polymer  coating  are  given  in 
equations  (59)  and  (60),  respectively;  they  are  repeated  here  as  equations  (A-l)  and  (A-2): 


r  ..  .  Tg-Tb 

dt  V c  ^Vc'kpA'  d 


(A-l) 


and 


It/  „  ^  d(Ta  +  Tb)  =  Tu  -  n  ..  v 

_  p  p  *  /Ip  [>  .  ~  AP  n(Tb  ~  Tx)« 

2  ot  a 

Equation  (A-l)  can  be  arranged  as 

(—  +  tv)  ra  -  Tc  T„  =  Qc, 
o  t 


where 


/r  p  Ac 

d  pc  CcVc 


Similarly,  rearranging  equation  (A-2)  yields 


Ta 


+  th  )  T„=  0, 


(A-2) 


(A-3) 


(A-4) 


(A-5) 


where 


2  Ap  _  2  h  A r 

d  pp  CP  VP  '  Th  Pp  CP  VP ' 


(A-6) 


Too  is  assumed  to  approach  zero,  and  d  denotes  the  thickness  of  the  coating.  These  constants  Ta, 
Tb,  Tc,  Th,  and  tp  have  units  of  inverse  time;  thus,  their  inverses  represent  a  characteristic  thermal 
relaxation  time  for  conduction  or  convection.  For  example,  l/rc  is  a  measure  of  the  relaxation 
time  for  conduction  of  heat  from  the  ceramic  to  the  coating. 

For  convenience,  write  the  operators  in  equations  (A-3)  and  (A-5)  as  a  coupled  system  of 
equations: 


A-l 


Dll  Ta  "  Dl2  Tb  -  QC1 


(A-7) 


D2\  Ta  +  D22  Tb  =  0, 

(A-8) 

where  the  operators  are 

b 

11 

+ 

0 

Di2  =  rc, 

°2' 

and 

D22=(%t  +  rp  +  Th). 

Multiplying  equation  (A-7)  by  D22and  multiplying  equation  (A-8)  by  D12 
combining  the  results  by  expressing  7*  in  terms  of  Ta  yields 

and  then 

(D,2  D„  +  D|2  D21  )Ta  =  D22  Qc. 

(A-9) 

Carrying  out  the  operations  in  this  equation,  using  the  definitions  of  the  D’s  and  r’s,  yields 

f  d2  d 

>r+s,&+£!  )T‘=c’e" 

(A- 10) 

where 

£\  =  2rc  +  rp  +  rh 

(All) 

£2  =(*’p  +  r*K-V|, 

(A- 12) 

*3  =*>  +**• 

(A- 13) 

Equation  (A- 10)  is  the  desired  energy  balance  equation;  it  is  equation  (61)  in  the  body  of  this 
report. 


A-2 


Electrical  Input  Power  (Watts) 


APPENDIX  B 

MEASURED  IN-AIR  DIELECTRIC  CONSTANT 
DISSIPATION  FACTOR  AND  POWER 


Figure  B-l.  Measured  Electrical  Input  Power  (W)  for  Cylinder  I, 
750  Vrms  (3.9  kV/in.)  at  5  kHz  in  Air  and  No  Coating 


Figure  B-2.  Measured  Dielectric  Constant  for  Cylinder  I, 
750  Vrms  (3. 9  kV/in.)  at  5  kHz  in  Air  and  No  Coating 


Figure  B-3.  Measured  Electrieal  Dissipation  Factor  for  Cylinder  /, 
750  Vrms  (3.9  kV/in.)  at  5  kHz  in  Air  and  No  Coating 


Figure  B-4.  Measured  Electrical  Input  Power  (W)  for  Cylinder  2, 
1000  Vrms  (5.4  kV/in.)  at  10  kHz  in  Air  and  No  Coating 


B-2 


Figure  B-7.  Measured  Electrical  Input  Power  (W)  for  Cylinder  J, 
1000  Vrms  (7.6  kV/in.)  at  10  kHz  in  Air  and  A 'o  Coating 


Figure  B-8.  Measured  Dielectric  Constant  for  Cylinder  3, 
1000  Vrms  (7.6  kV/in.)  at  10  kHz  in  Air  and  No  Coating 
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